Search for: All records

Abstract Disk continuum reverberation mapping is one of the primary ways we learn about active galactic nuclei (AGN) accretion disks. Reverberation mapping assumes that time-varying X-rays incident on the accretion disk drive variability in UV–optical light curves emitted by AGN disks and uses lags between X-ray and UV–optical variability on the light-crossing timescale to measure the radial temperature profile and extent of AGN disks. However, recent reverberation mapping campaigns have revealed oddities in some sources, such as weakly correlated X-ray and UV light curves, longer than anticipated lags, and evidence of intrinsic variability from disk fluctuations. To understand how X-ray reverberation works with realistic accretion disk structures, we perform 3D multifrequency radiation magnetohydrodynamic simulations of X-ray reprocessing by the UV-emitting region of an AGN disk using sophisticated opacity models that include line opacities for both the X-ray and UV radiation. We find there are two important factors that determine whether X-ray irradiation and UV emission will be well-correlated: the ratio of X-ray to UV luminosity and significant absorption. When these factors are met, the reprocessing of X-rays into UV is nearly instantaneous, as is often assumed, although linear reprocessing models are insufficient to fully capture X-ray reprocessing in our simulations. Nevertheless, we can still easily recover mock lags in our light curves using software that assumes linear reprocessing. Finally, the X-rays in our simulation heat the disk, increasing temperatures by a factor of 2–5 in the optically thin region, which could help explain the discrepancy between measured and anticipated lags. 

Abstract Active galactic nuclei (AGN) light curves observed with different wave bands show that the variability in longer wavelength bands lags the variability in shorter wavelength bands. Measuring these lags, or reverberation mapping, is used to measure the radial temperature profile and extent of AGN disks, typically with a reprocessing model that assumes X-rays are the main driver of the variability in other wavelength bands. To demonstrate how this reprocessing works with realistic accretion disk structures, we use 3D local shearing box multifrequency radiation magnetohydrodynamic simulations to model the UV-emitting region of an AGN disk, which is unstable to the magnetorotational instability and convection. At the same time, we inject hard X-rays (>1 keV) into the simulation box to study the effects of X-ray irradiation on the local properties of the turbulence and the resulting variability of the emitted UV light curve. We find that disk turbulence is sufficient to drive intrinsic variability in emitted UV light curves and that a damped random walk model is a good fit to this UV light curve for timescales >5 days. Meanwhile, X-ray irradiation has negligible impact on the power spectrum of the emitted UV light curve. Furthermore, the injected X-ray and emitted UV light curves are only correlated if there is X-ray variability on timescales >1 day, in which case we find a correlation coefficientr= 0.34. These results suggest that if the opacity for hard X-rays is scattering dominated as in the standard disk model, hard X-rays are not the main driver of reverberation signals. 

Abstract The variability of quasar light curves can be used to study the structure of quasar accretion disks. For example, continuum reverberation mapping uses delays between variability in short and long wavelength bands (shortlags) to measure the radial extent and temperature profile of the disk. Recently, a potential reverse lag, where variations in shorter wavelength bands lag the longer wavelength bands at the much longer viscous timescale, was detected for Fairall 9. Inspired by this detection, we derive a timescale for theselongnegative lags from fluctuation propagation models and recent simulations. We use this timescale to forecast our ability to detect long lags using the Vera Rubin Legacy Survey of Space and Time (LSST). After exploring several methods, including the interpolated cross-correlation function, a Von-Neumann estimator,javelin, and a maximum-likelihood Fourier method, we find that our two main methods,javelinand the maximum-likelihood method, can together detect long lags of up to several hundred days in mock LSST light curves. Our methods work best on proposed LSST cadences with long season lengths, but can also work for the current baseline LSST cadence, especially if we add observations from other optical telescopes during seasonal gaps. We find that LSST has the potential to detect dozens to hundreds of additional long lags. Detecting these long lags can teach us about the vertical structure of quasar disks and how it scales with different quasar properties. 

ABSTRACT Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination i < 90°) decrease inclination, as well as semimajor axis (a) and eccentricity (e) until orbital alignment with the gas disc (‘disc capture’). Captured stellar-origin black holes (sBH) add to the embedded AGN population that drives sBH–sBH mergers detectable in gravitational waves using LIGO–Virgo–KAGRA or sBH–SMBH mergers detectable with Laser Interferometer Space Antenna. Captured stars can be tidally disrupted by sBH or the SMBH or rapidly grow into massive ‘immortal’ stars. Here, we investigate the behaviour of polar and retrograde orbiters (i ≥ 90°) interacting with the disc. We show that retrograde stars are captured faster than prograde stars, flip to prograde orientation (i < 90°) during capture, and decrease a dramatically towards the SMBH. For sBH, we find a critical angle iret ∼ 113°, below which retrograde sBH decay towards embedded prograde orbits (i → 0°), while for io > iret sBH decay towards embedded retrograde orbits (i → 180°). sBH near polar orbits (i ∼ 90°) and stars on nearly embedded retrograde orbits (i ∼ 180°) show the greatest decreases in a. Whether a star is captured by the disc within an AGN lifetime depends primarily on disc density, and secondarily on stellar type and initial a. For sBH, disc capture time is longest for polar orbits, low-mass sBH, and lower density discs. Larger mass sBH should typically spend more time in AGN discs, with implications for the spin distribution of embedded sBH. 

Abstract While most simulations of the epoch of reionization have focused on single-stellar populations in star-forming dwarf galaxies, products of binary evolution are expected to significantly contribute to emissions of hydrogen-ionizing photons. Among these products are stripped stars (or helium stars), which have their envelopes stripped from interactions with binary companions, leaving an exposed helium core. Previous work has suggested these stripped stars can dominate the Lyman Continuum (LyC) photon output of high-redshift, low-luminosity galaxies post-starburst. Other sources of hard radiation in the early universe include zero-metallicity Population iii stars, which may have similar spectral energy distribution (SED) properties to galaxies with radiation dominated by stripped-star emissions. Here, we use four metrics (the power-law exponent over wavelength intervals 240–500 Å, 600–900 Å, and 1200–2000 Å, and the ratio of total luminosity in FUV wavelengths to LyC wavelengths) to compare the SEDs of simulated galaxies with only single-stellar evolution, galaxies containing stripped stars, and galaxies containing Population iii stars, with four different initial mass functions (IMFs). We find that stripped stars significantly alter SEDs in the LyC range of galaxies at the epoch of reionization. SEDs in galaxies with stripped stars have lower power-law indices in the LyC range and lower FUV to LyC luminosity ratios. These differences in SEDs are present at all considered luminosities ( M UV > − 15 , AB system), and are most pronounced for lower-luminosity galaxies. Intrinsic SEDs as well as those with interstellar medium absorption of galaxies with stripped stars and Population iii stars are found to be distinct for all tested Population iii IMFs. 

Abstract High-resolution numerical simulations including feedback and aimed at calculating the escape fraction (f_esc) of hydrogen-ionizing photons often assume stellar radiation based on single-stellar population synthesis models. However, strong evidence suggests the binary fraction of massive stars is ≳70%. Moreover, simulations so far have yielded values off_escfalling only on the lower end of the ∼10%–20% range, the amount presumed necessary to reionize the universe. Analyzing a high-resolution (4 pc) cosmological radiation-hydrodynamic simulation, we study howf_escchanges when we include two different products of binary stellar evolution—stars stripped of their hydrogen envelopes and massive blue stragglers. Both produce significant amounts of ionizing photons 10–200 Myr after each starburst. We find the relative importance of these photons to be amplified with respect to escaped ionizing photons, because peaks in star formation rates (SFRs) andf_escare often out of phase by this 10–200 Myr. Additionally, low-mass, bursty galaxies emit Lyman continuum radiation primarily from binary products when SFRs are low. Observations of these galaxies by the James Webb Space Telescope could provide crucial information on the evolution of binary stars as a function of redshift. Overall, including stripped stars and massive blue stragglers increases our photon-weighted mean escape fraction () by ∼13% and ∼10%, respectively, resulting in. Our results emphasize that using updated stellar population synthesis models with binary stellar evolution provides a more sound physical basis for stellar reionization. 

Abstract The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA’s first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultra-compact stellar-mass binaries, massive black hole binaries, and extreme or interme-diate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe.